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1.  Introduction 


We  eouidor  an  infinite  aorror  queue  to  which  cuatoaora  arrive  according 
to  a  non-homogeneous  Poiaaon  Procaaa  CHHPP)  with  known  intensity  A(t) ,  t  >_  o. 

A  customer  arriving  at  tine  a  innadlately  enters  service  and  with  probability 
F(s,w)  will  conplete  service  by  tine  w  where,  of  course,  F(s,w)  -  o  for 
w  <  a.  Service  tinea  for  custenera  era  aasuned  to  be  statistically  indepen¬ 
dent  of  each  other.  Such  a  queue  will  be  referred  to  as  an  M(t)/G(t)/** 
queue  to  indicate  that  the  arrival  process  is  a  HHPP,  that  a  customer's  service 
tine  nay  depend  on  the  tine  ha  entered  service  and  that  there  are  an  infinite 
nunber  of  servers.  If  F(s,w)  ■  G(v-s)  is  a  function  only  of  the  difference 
v-a  and  A(t)  -  A  for  all  t,  we  have,  aa  a  special  case,  the  familiar  M/G/« 
queue. 

An  Important  stochastic  process  of  interest  in  queueing  theory  is  the 
counting  process  (D(t),  t  >.  o>  where  D(t)  is  the  nunber  of  departures  from 
the  queue  by  tine  t.  In  thia  note  we  establish  and  use  a  splitting  theorem 
for  HHPP  to  show  that  for  the  M(t)/C(t)/«*  queue,  {D(t),  t  >  o}  is  a  NHPP. 
Furthermore,  we  show  that  the  norther  of  departures  from  the  queue  during  any 
finite  tine  interval  is  independent  of  the  nusher  of  customers  still  in  service 
at  the  end  of  that  Interval.  As  special  cases  of  this  result  we  show  that  the 
transient  output  of  an  M/G/**  queue  is  a  HHPP  and  that  In  steady  state  the 
output  process  of  an  M/G/**  queue  is  a  stationary  Poisson  Process  with  rate 
equal  to  the  rate  of  the  Poisson  input  process.  These  results  for  the  M/G/** 
queue  were  first  proved  more  laboriously  by  Mlrasol  [4] . 

Hlllestad  and  Carillo  [1]  have  shown  that  the  nunber  of  busy  servers  at 
tine  w  in  the  M(t)/G(t)/**  queue  is  Poisson  distributed  with  naan  depending 
on  w  and  the  form  of  the  service  time  distribution.  This  result  can  be  simply 
obtained  by  judiciously  using  the  splitting  theorem  established  in  Section  3. 

To  the  author's  knowledge,  no  study  has  been  nsde  of  the  output  process 
of  the  M(t)/G(t)/**  or  M(t)/G(t)/s  queues.  In  a  sequel  to  this  report,  we 
shall  study  the  output  proesss  of  M(t)/M/s  queues. 

2.  Mon-Homogeneous  Poisson  Processes 

Definition  2.1.  Let  {H(t),t  >  o)  be  a  stochastic  counting  process  and 
let 

a.  H(o)  ■  o 

b.  {N(t),  t  >  o>  have  Independent  lncrsnents. 
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c.  ?r(2  or  more  events  in  (t,t+h))  -  0(h). 

d.  Pr (exactly  one  event  in  (t,t+h))  »  ^ (t)h  +  0(h). 

Then{N(t),  t  >_ o)is  said  to  be  e  non-hoaogeneous  Poisson  Process  with 

Intensity  function  X(t). 

Let  M(t)  -  /x(e)ds.  Then  it  is  straight forward  to  show  that  Definition 
2.1  is  equivalent  to  Definition  2.2. 

Definition  2.2.  (H(t),  t  >_  o)is  said  to  be  a  HHPP  with  mean  value  func¬ 

tion  (KVF)  M(t)  if 
a'.  N(o)  ■  o 

b'.  (H(t) ,  t  ^  o}  has  independent  increments, 

c* .  for  all  s  <  t  and  n  £  o 

Pr(N(t)  -  H(s)  -  n)  -  «“CM(t)  -  H(s))  .(M(t},_-_M(s))_ 

In  this  report,  we  shall  assume  the  HVF  of  all  HHPP  to  be  completely 
differentiable.  The  proofe  presented  here  can  be  extended  to  the  case  where 
the  MVF  is  not  completely  differentiable  but  this  would  add  unnecessary 
technical  complexity  to  the  discussion  herein. 

3.  A  Splitting  Theorem 

It  is  well  known  (see  IS]  or  16])  that  if  {N(t),  t  ^  o}  is  a  stationary 
Poisson  Process  with  rate  X  and  an  event  that  occurs  at  time  t  is  classified 
as  type  1  with  probability  P^,  1*1, . . .n,  independently  of  the  classification 
of  other  events  (ZP^  ■  1),  than  the  stochastic  processes  {H^(t),  t  ^  o},  1*1,..., 
n  ere  Independent  stationary  Poisson  Processes  with  ratal  xP^,  1  *l,...,n.  In 
this  section  we  prove  an  analogous  result  for  HHPP. 

Theorem  1:  Let  (H(t),  t  >  o)  b«  a  HHPP  with  intensity  X(t)  and  MVF  M(t). 

An  event  that  occurs  at  time  w  >  e  li  classified  as  a  type  1  event  with 
probability  P1(w)  and  as  a  type  2  event  with  probability  ?2(w)  -  1  -  P^(w). 
Classification  of  each  event  is  independent  of  the  classification  of  other 
events.  Let  H^(t)  be  the  nuaber  of  type  1  events  that  have  occurred  by 
time  t.  Then  {H^(t),  t  >  o}  is  a  HHPP  with  intensity  X(t)P^(t)  and  MVF 
Mj^t)  *  f  X(s)Pi(s)ds»  i  -  1,2.  Furthermore,  {^(t),  t  >_  o)  and  {H2(t),t  >_  o} 
are  statistically  Independent. 

Before  proceeding  with  the  proof  of  Theorem  1,  we  need  to  establish  the 
following  Leans  which  is  given  as  an  exercise  in  both  Boss  [6,p29]  and 
Person  [5,pl43]. 
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Laana  1: _ Let  {H(t),  t  *,  o>  be  a  HHPP  with  Intensity  X(t)  and  MFF  M(t). 

Given  that  H(t)  -  a,  tha  Joint  diatrlhatloa  of  tha  a  apocha  at  which  tha  events 

occurred,  x^  ,  t2 . t&  <  t  li  tha  aane  aa  if  they  ware  tha  order  statistics 

corresponding  to  n  Independent  raadon  variables  y^,  y2,...yn  with  canon 
distribution  function 

F<«>  "  Iftf  0  <  s  <  t  (1) 

Proofs  Lot  <,  y j2]  <. .  .<  y^j  bo  the  order  statistics  corresponding 

to  7i»  7 2’ Then  clearly,  the  joint  density  of  these  order  statistics 
Is  given  by 

g -  n!f(v1)f(v2)...f(vn),v1  1  »2  £  ...<, 

Where 

«(»)  -  &  «*>  - 

Thus,  n 

*r<ym  -  w.,  i  -  1,2, . .  .n>— *  X(w.)  (2) 

I1]  1  K(t)n  1  1 

How  we  show  that  the  joint  distribution  of  the  ...,1^  given  H(t)  -  n.  Is 
identical  to  (2) .  Let  iv2  £•••  w#  <.  t  and  leth^  be  snail  enough  so  that 
w^  +  h^  <  wi+2.  Then 

Pr(wt  <t1<*1  +  h2,  1-1,2, ...,n|n(t)  -  n) 

-  Pr(wt  <  <  wt  +  ht,  i  -  1,2,..., n)/Pr(H(t)  -  n) 

-  Pr  (  HCWj  +  h2)  -  I<w4)  -  1,  1-1,..., n ; 

H(w^)  -  *Cw1-1-Hi1__1)  -  o,  1-2,..., np 

■(Wj)  -  H(o>  »  o; 

H(t)  -  N(w  +  h  )  -  o)/Pr(H(t)  -  n)  (3) 
n  n 

Frost  c'  of  Definition  2.2, 

Pr(JI(w1)  -  N(w±_1  +  ht-1)  -  o)  -  e-0*^  "  M(wi-1  +  Vl^ 

PrOKwj)  -  H(o)  -  o)  -  e-*^ 

Pr(H(w1  +  h±)  -  H(w1)  -  1)  -(M(w14h1)-M(w1))  .-OKVV  -  M(w1)) 

PrQl(t)  -  H(wn+hft)  -  o)  -  e"Wt)  “  (4) 
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Sines  the  interval*  Iv^v^+hj]  <lo  not  overlap,  ve  have  from  (b)  of  Definition 
2.2  and  (4)  that  (3)  bee one*  after  sons  reduction 

#-M(t)  J  (M(v1+h1)  -  M(vt)) 


-M(t)  H(t)n 
*  n! 


•  nl _  n 

M(t)n  l  O^VV  "  M<»i» 


Therefore, 


Pr^wl£cl<*i4fci»  1-1 »  »|*(t)-n)  n,  n  M^-fh^  -  M(w^) 


(5) 


v  h. 


M(t)®  1 


HOw  take  the  linit  of  both  eidee  as  the  h^  ♦  o,  i*l. . .,n  uniformly.  Since 


11a  M(w^*Hi^)  -  M(Wj) 


h^-K) 


MCw1) 


X(w±) 


and  since  the  HVPs  are  differentiable,  the  limit  of  the  right  hand  side  of  (5) 
exists  which  therefore  Implies  the  limit  on  the  left  side  of  equation  (5)  also 
exists  and  it  is,  in  fact,  an  ordinary  probability  density  function.  Thus, 

fGrj . wtt-|N(t)  -  n)  -  lia  PrCw^^^^.i-l. ..,n|K(t)-n) 


hj+o 

1*1,..., n 


$ 

1 


;  *<*i> 


which  is  precisely  the  probability  density  function  (2). 

Leans  1  states  that  given  N(t)  ■  n,  the  n  event  times,  ^i»^2*****Tn  w*'*n 
considered  as  unordered  randon  variables  have  the  same  joint  distribution  as 
n  Independent  random  variables  with  common  distribution  F(s)  given  by  (1) . 
Therefore,  so  a  consequence  of  Lemma  1,  given  that  we  know  an  event  occurred 
in  (o,t),  the  actual  time  that  the  event  occurred  has  distribution  function  F(s) 
sad  is  Independent  of  the  time  of  the  other  events  occurring  in  (o,t). 

5 


Proof  of  Theoraa  1.  VO  prove  Theoraa  1  bp  first  showing  that  (N±(t) ,  t  >_  o> 
i  ■  1,2  satisfy  Definition  2.2  for  RRPP.  We  shell  show  that  F^t)  Is  Poisson 
distributed  with  assn  /i(s)p^(s)ds  by  eslculsting  the  generating  function  of 
N^(t) ,  1  ■  1,2.  The  independence  of  the  two  processes  Is  established  by  cal¬ 
culating  the  joint  generating  function  of  H^(t)  and  N2(t)  and  showing  that 
this  joint  generating  function  Is  the  product  of  the  Individual  generating 
functions. 

Clearly,  N^(o)  ■  NgCo)  ■  o  so  that  (a')  of  Definition  2.2  holds.  For 
1  -  1,2,  t  >_  o  and  |*|<1,  let 

VO 

g1(z,t)  -  E[*  1  ] 

be  the  generating  function  for  N^(t) .  By  conditioning  on  N(t)  and  applying 
the  Law  of  Total  Probability  we  have  that 

VO 

8l(*,t)  -  S^jlEC*  X  |  I(t) (6) 

Proa  Leaaa  1,  given  that  an  event  occurred  in  (o,t],  the  probability  It 
occurred  at  ties  z,  o  <  x  <  t  Is  .  Given  an  event  occurred  at  tine  x, 

the  probability  It  was  classified  as  a  type  1  event  is  p1(x) .  Therefore, 
given  that  an  event  occurred  In  (o,t]  the  probability  It  was  a  type  1  event 


Independently  of  the  other  events  that  occurred  In  (o,t].  Therefore,  given 
H(t),  VO  la  a  binoalal  randoa  variable  with  paraaeters  N(t)  and  y±  so  that 
the  conditional  generating  function  of  R^(t)  le 

E[*  1  |N(t)]  -  (1-y±  +  Yl*)H(t) 

and  using  this  in  (6)  we  have 

«iU.O  -  K„(t)  ICi  -  y±  + 

Since  R(t)  la  Polaaon  distributed  with  assn  M(t)  and  1 1— I  < 1  • 

-  .-««  (Ti-ITi*)  .  .-M<t)Yl(l-.)  (8) 


< 


But  (8)  la  Just  eha  generating  function  of  a  Polaaon  random  variable  with  aean 


YjlKt)  K*)  P1(x)dx  (9) 

O 

Thus,  N^(t)  1  »  1,2  la  Polaaon  distributed  with  naan  given  by  (9)  for  all 
t  >.o.  By  a  similar  argument  one  can  show  that  fgr  any  Interval  (t^.tj], 

(tj)  -  N^(t^)  la  Polaaon  dlatrlbuted  with  moan  /2X(z)p^(x)dz  since  given 

that  an  event  occurred  in  (t^,t2J»  It  la  type  1  with  probability 

\  Ctj^.tj)  -  J2  p±(x)  -  M(t-))  1-1.2 

*1 

Independently  of  the  other  eventa  occurring  In  (t^,^]. 

That  each  proceaa  baa  Independent  lncramenta  follow  straightforwardly 
from  the  fact  that  (H(t)  •  t  >.  o>  ha a  Independent  lncramenta  and  that  each 
event  la  typed  Independently  of  other  eventa.  If  we  consider  any  finite  number 
of  non-over lapp lng  time  Intervale,  the  number  of  type  1  events  that  occur  In 
each  Interval  depends  only  on  the  nuaber  of  events  of  the  process  (N(t),  t  >_  o} 
that  occurred  In  that  Interval  and  tha  typing  mechanism  since  only  events  of 
the  process  (N(t) ,  t  >  o}  occurring  la  a  particular  Interval  can  be  classified 
and  counted  as  a  type  i  event  in  that  interval.  Since  {N(t) ,  t  ^  o}  has 
Independent  increments,  the  fact  that  the  processes  {N^t),  t  >_  o}  1-1,2 
have  Independent  increments  follows  immediately  by  conditioning  on  the  nuaber 
of  events  of  (N(t),  t  _>  o}  that  occurred  in  each  Interval. 

All  the  conditions  of  Definition  2.2  have  been  satisfied  and  therefore 
{N^(t),  t  >_  o}  1-1,2  are  NHPP. 

All  that  remains  now  In  proving  Theorem  1  Is  to  show  that  ^(t)  and 
H2(t)  are  Independent  for  all  t.  Consider  the  joint  generating  function 

gCaj.^.t)  -  E(*1Ill(t)*2S2(t)l 

m  m 

-  S  I  {« J  PrOi^t)  -  n,  H2(t)  -  m|  H(t)  -  n  +  m) 

n— o  w— o 

•  Pr(N(t)  -  m  +  m)> 

From  our  previous  discuss Iona  we  realise  that 

PrCHjtt)  -  n,  H2(t)-  *|H(t)  -  *  +  n)  ^ 


since  given  N(t)  ■  a  +  n,  there  ere  (  )  veye  of  choosing  n  events  to  be 

type  1  (end  thereby  choosing  the  *  to  be  type  2) .  Thus, 

g<«rVt)  -  l  O  Y*  e-M(t>M(t)n+“/(n^) 

n«o  * 


.  Z  M(t)n  l  «-M<t)lf2  H(t)» 


n«o 


n! 


m! 


where  we  have  used  the  fact  that 

e-M(t)  .  e-M(t)Y1  e-M(t)Y2  Y1  +  Y2  "  1 

Thus, 

g(*ltz2tt)  -  gOtj^t)  g(s2,t) 

Since  this  holds  for  all  z^,  s2  and  t,  N^(t)  and  N^Ct)  are  independent  and 
this  completes  the  proof. 

Note  that  It  is  trivial  to  extend  Theorem  1  to  the  case  when  an  event 
at  time  t  can  be  classified,  independently  of  other  events,  into  one  of  n 
categories  with  the  probability  that  the  event  Is  classified  as  type  1  being 
p^Ct)  with  E  p^Ct)  *  1  for  all  t.  This  splitting  theorem  for  NHPP  Is  useful 
In  proving  many  Interesting  properties  of  NHPP.  In  particular,  we  shall  use 
It  to  show  that  the  output  process  of  the  M(t)/G(t)/aD  queue  is  Indeed  a  NHPP. 

4.  The  Output  of  an  M(t)/0(t)/"  Queue 

In  this  section  we  will  determine  the  distribution  of  the  number  of  depar¬ 
tures  In  an  arbitrary  Interval  (t,  t+T)  of  length  T  In  the  M(t)/G(t)/">  queue. 
Recall  that  In  this  queue,  service  times  are  Independent  of  each  other  and 
that  If  a  customer  arrives  at  time  s,  the  probability  he  will  complete  service 

by  time  w  la  7(s,w),  with  /  d  (F(s,w))  -  1. 

s  w 

Let  F(s,w)  -  1  -  F(s,w) .  We  shall  classify  any  arrival  Into  one  of  three 
types: 

Type  1:  if  the  arrival  will  depart  the  system  In  (t,t+T] 

Type  2:  If  the  arrival  will  depart  the  system  In  (t+T,*) 

Type  3:  If  the  arrival  will  depart  the  system  In  (o,t]. 


Let  p^Cx)  be  the  probability  that  c  customer  arriving  at  tine  x  la  classified 
as  a  Type  1  customer*  1  *  1*2*3.  Then*  for  all  x  >_  o 

Pj^Cx)  -  F(x,t+T)  -  F(x»t) 

P2(x)  -  F(x,t?*-T)  -  1  -  F (x* t+T) 

P3(x)  -  F(x,t) 

where ^we  have*  of  course*  adopted  the  convention  that  F(x»y)  •  o  if  x  >  y.  Rote 
that  2  j.  > (x)  -  1  for  all  x  and  since*  by  assumption*  each  arrival  Is  classified 
Independently  of  other  arrivals*  the  conditions  of  Theorsm  1  have  been  satisfied. 
Therefore*  If  N^(w)  -  the  number  of  customers  mho  have  arrived  by  time  w  and 
are  categorized  as  Type  1*  we^have  from  Theorsm  1  that  (H^Cv)*  w  ,>  o)  , 

1  -  1*2*3  are  NHPP  with  MVF  A(s)p1(s)ds,  1  -  1*2*3*  and  ere  mutually 
Independent.  We  recognize  tfiat  {N^OO ,  w  >.  o)  is  simply  the  counting  process 
describing  departures  from  the  queue  In  (t,t+T] .  In  particular*  Ij(tfl)  Is 
the  number  of  arrivals  by  t+T  who  have  departed  in  (t,t+T]  end  N^(t+T)  Is 
Independent  of  N2 (t+T)  which  Is  the  number  of  arrivals  by  t+T  that  are  still 
in  the  system  at  time  t+T.  Thus*  he  have  the  Important  and  Interesting  result 
that  for  an  M(t)/G(t)/<*>  the  number  of  departures  la  any  Interval  la  Independent 
of  the  number  remaining  In  the  system  at  the  and  of  the  Interval.  This 
generalizes  Mlrasol's  result  for  H/G/»  queues.  Furthermore*  H^(t+T)  la  Inde¬ 
pendent  of  Nj(t+T)  -  N3 (t)  which  la  the  number  of  departures  In  (o*t] .  In 
fact*  If  we  take  any  countable  number  of  non-overlapping  Intervals*  using  the 
same  analysis  as  above*  It  Is  straightforward  to  sheer  that  the  departures  during 
any  of  these  Intervale  la  Poisson  distributed  and  la  Independent  of  departures 
during  any  other  Interval.  (If  there  are  R  non  over-lapping  Intervals  let 
p  (x)  be  the  probability  that  an  arrival  at  time  x  will  complete  service  In 

®  .L 

the  n  <_ N  interval.  Then  argue  as  above).  Clearly  than*  the  conditions 
of  Definition  2.2  are  satisfied*  and  by  letting  t*o  we  see  H^(T)  -  D(T)  so 
that  we  have  established  the  fact  that  (D(w)  *  w  >.  o)  la  a  NHPP  and  for  all 
T  >.  o*  D(T)  Is  independent  of  the  nuhber  still  In  service  at  time  T.  Also* 
we  have  as  a  by-product  that  N2(t+T)  la  Poisson  distributed  so  that  the  number 
of  busy  servers  at  time  t  in  the  M(t)/G(t)/*»  queue  Is  Poisson  distributed 
with  mean  /  *(,)  F(*,t)dx. 

For  ail  tio  then*  D(t)  is  Poisson  distributed  with  MVF 

t 

(t)  -  /  *(*)  F(x*t)dx 
o  9 


(10) 


and  intensity  function 

V°  "  ft  *D(t)  •  /  *<*>  *<*>*>**  +  x(t)  F(t.t)  (U) 

r.  o 

ihKi  f  (x»t)  •  and  ?(t,t)  -  o  only  If  the  service  distribution  at 

tin  t  has  an  atom  at  o. 

Betumiag  to  the  process  {K^ (v) ,  v  >_  o}  ve  note  that  for  all  t,T  >_  o 


Hhjft.T)) 


.Quaber  of  arrivals  by  t+T  who  . 

1  departed  the  systen  in  (t,t+Tr  . 


■  E  [number  of  departures  in  (t,t+T]] 
t+T  t+T 

-  /  X(x)p,(x)dx  •  /  X(x)[F(x,t+T)  -  P(x,t)]dx 

o  o 

t 

-  /  X(x)  [f Cx»t+T)  -  F(x,t)]  dx 

°  t+T 

+*  /  X(x)  F(x,t+T)dx  (12) 

t 

The  first  tern  in  (12)  can  easily  be  seen  to  be  the  number  of  arrivals  in 
(o,t]  who  depart  in  (t, t+T]  and  the  second  tern  la  the  number  of  arrivals  in 
(t,t+T]  uho  igg|| rt  in  (t,t+T] .  (12)  can  also  be  obtained  by  noting  that 

Efo^t+T)]  -  /  XJ)(x)dx. 

5.  Special  Cases 

a.  M(t)/G/«»  Queue 

For  w  >.  s  let  F(s,v)  ■  G(vs)  be  a  function  only  of  the  difference  vs 
This  la  just  the  classical--  definition  of  service  tinea  in  queues.  From  (12) 
(or  equivalently,  from  (11))  ve  have  that  H1(t+T)  la  Poisson  distributed  with 

mean 

t  t+T 

/X(x) [G(t+T-x)  -  G(t-x)]dx  +•  /  x(x)  G(t+T-x)dx 
o  t 


„  t+T  t 

/  X(x)  G(t+T-x)dx  -  /X(x)G(t-x)dx  (13) 

o  o 

b.  H/8/s  Queue 


Let  the  service  tine  distribution  be  as  in  Case  a*  and  let  X(t)  -  X 


for  all  t>_  o.  From  (13)  we  see  that 

t+T 

E[N,(t4T)]  -  /  XG(y)dy  (14) 

t 

From  (14)  we  note  that 

t+T 

11m  E[N-(t+T)]  »  11m  /  XG(y)dy  -  XT 

t-+«  t 

so  that  in  steady  state,  the  expected  number  of  departures  over  an  Interval 
of  length  T  Is  simply  XT.  Thus,  we  would  expect  that  for  the  M/G/«  queue, 
the  steady  state  output  rate  would  equal  the  Input  rate  X .  This  Is  verified 
directly  from  (11)  since  the  Intensity  function  of  the  departure  process  is 


Xj^t)  ■  /  Xg(t-x)dx  +  XG(o) 


-  XG(t) 


lim  ^_(t)  -  x 

t-+«> 


This  interesting  result  for  M/G/«  queues  was  first  obtained  by  Mirasol  [5] . 


Let  F(x 


queues 

(o  if  t  <  X  + 
.t)  -  < 

(1  O.W. 


for  all  x  represent  the  distribution  for  a  deterministic  service  time  of  length  T. 

t  . 

From  (10),  11. (t)  -  /  X(x)dx  -  the  expected  number  of  arrivals  in  (t-T,t) 

"  t-T 


and  from  (11) 


Xjj(t)  -  X(t-T) 


That  is,  the  departure  intensity  at  t  is  simply  the  arrival  intensity  T  time 
units  previous.  This  is  as  it  should  be  since  for  h  small, 

Fr(daparture  in  (t,t-H0) 

•  Fr (arrival  in  (t-T,  t-T+h) 

-  X(t-T)h  +  0(h) 

by  our  assumptions  on  the  arrival  process. 


▲  HHEP  with  Intensity  function  X(t)  la  said  to  bo  "dying"  If 

(1)  o  i  X(t)  <  •  for  all  t  >.  o 

(2)  11b  X(t)  -  o 

Thus,  a  NHPP  la  dying  if  its  iatanalty  function  la  dying  out  with  tlae  and  if 
during  any  finlta  Interval  the  expected  nuBber  of  events  of  that  process  that 
occur  la  finite.  Intuitively*  we  would  expect  that  if  the  input  process  to  a 
queue  la  dying  that  the  output  proceae  should  be  dying  as  wall.  Theorem  2, 
a  Tauber lan  type  theor«a,  verifies  our  intuition. 

Theorem  2:  For  the  M(t)/Q(t)/"  queue*  the  departure  process  is  dying 
if  and  only  if  the  arrival  process  is  dying. 

Proof:  Ws  shall  prove  Theorem  2  for  queuee  whoee  arrival  intensity 
function  and  its  derivatives  are  Lsplsce  Transformable  for  all  la(a)  1  o. 

(See  [7]).  let  \(*)t  *D(t)  be  the  intensity  functions  of  the  errlval  and 
departure  processes*  respectively,  from  (11) 

t 

11m  X  (t)  ■  11m  /  X  (n)f  (x*t)dx  +  11m  X  (t)F(t,t) 

f**  (4*  O  *  t**  * 

t 

■  /  lA(x)f(x,t)dx  (15) 

O 

since  o  <  11m  X.(t)F(t*t)  <_  11m  X. (t)  ■  o.  Therefore*  to  study  the  limiting 
behaviour  of  ^(t)  it  sufflceito  study  the  limiting  behaviour  of  the  function 

t 

b(t)  -  /  X  (m)f(x*t)dx  (16) 

o 

let  XA(a)  -  /  e“**XA  Ct)dt_ 

*  0  • 

X0(e)  -  /  •"•\Ct)dt 
o 

f(x,s)  -/."•  (t“^)f(x.t)dt 
X 

h(s)  -  r i"*1  b(t)dt 


Taking  tha  Laplaca  Tranafor*  (LT)  of  both  aldaa  of  (16)  wa  aaa 


Lot 


b(a) 


/  a 

o 

• 

/ 
o 


-at  | 


AA(*)f(*,t)dat|  dt 

•“•<t"*)f(*.t)dtj  dx 

*  -«  * 

■  /  •  XA(x)f<x,a)dx 


(17) 


Y  (•)  ■  *«P  f(x,a) 


Y  (a)  ■  Inf  f (x,a) 


Bacall  that  for  find  x,  /  f(x,y)dy  -  1  a  Inca  It  rapraaanta  tha  (propar)  aarvice 
tine  distribution  of  a  cultowar  antarlng  aarrlca  at  tlaa  x.  Thar afore,  for 
all  o  <a  <• 

o  <  X  (a)  <  y  (a)  <,  1 
Fro*  (17)  than  wa  hava  that 

X  (■)  V-)  i  b(a)  <  XA(a)  y  (•) 

and  tbarafora, 

•  X  (■)  *A(«)  1 •  b  (•)  <  li*+  a  y  (•)  XA(a)  (18) 
a4©  r*o  o-*o 

Front  (18)  wa  can  anally  prove  Thaoraa  2  alnca  by  tha  Final  Value  Thaoran  for 
LT  wa  know 

^A(t)  -  li*+  a  lA(a) 

0*0 

b(t)  -  11*.  a  b  (a) 

9+0 

Wow,  from  (18)  wa  have  that 

11*.  ax. (a)  <  11*.  a  b  (a)  <  11*.  ax.  (a)  lia.TCa)  (1! 

040  9+0  ~  040  *  040 
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11* 

fc4»» 

and 

11* 

t+m 

Say  11* 

tr*» 

XA(t)  -  o. 

11*+  ^(a) 

and  sines  11*.  £  (■)  ■  11*.  y  (■)  ■  1  ns  bars,  by  applying  the  Final  Value 
Theorem,  **° 

o  £  11a  b(t)  £  o 

”  fri¬ 
end  thus  11a  b(t)  ■  11a  l_(t)  ■  o 

fM»  f***  ° 

Fro*  the  left  Inequality  la  (19)  vs  also  sea  that  If  11a  XD(t)  ■  o  it  wist  be  that 
11*+  s  XA(s)  -  o  and  ths  theoraa  Is  established.  tr*m 

®"*°  The  techniques  developed  In  the  proof  of  Theoraa  2  lead  to  an  Interesting 
result  for  gensral  Lap lacs  transformable  functions  that  can  be  represented 
as  the  convolution  of  two  functions. 

Theorem  3:  Let  b(t)  be  Laplace  Transformable  and  let 

t 

b(t)  •  /  g(x)h(t-x)dx  (20) 

o 

be  the  convolution  gf  two  non-negative,  non- identically  aero,  proper  functions 
g  and  h  such  that  /  g(x)dx  »  •  •  /  h(x)dx.  Then  11a  b(t)  ■  o  if  and  only 
if  .  °  °  tr*~ 

11a  g(a)  -  11*  h(s)  -  o 

Proof:  Since  b(t)  Is  a  convolution, 
b(a)  -  i(s)h(s) 

AAA 

where  b(s),  g(s),  h(a)  are  the  LT  of  b,  g,  and  h,  respectively.  Then  by  the 
Final  Value  Theorem  we  have 

11a.  sb(s)  «  (11a,  sg(s))(llaf  h(s))  -  (11a  g(s))(lla.  sh(s)) 
o*o  w*o  w*o  s-*o  e*o 

»  •  A  • 

Since  11a.  g(s)  ■  /  gOOdx  ■  •  •  11a.  h(s)  -  /  h(x)dx 

e*o  o  o*oT  o 

Theoraa  3  follows  laasdlately. 

Inspection  of  the  convolution  Integral  (20)  Intuitively  shows  why  Theorem  3 
Is  true.  If  both  functions  g(x)  and  h(x)  did  not  go  to  sero  sows  significantly 
positive  mass  would  be  accuanlated  when  Integrating  fro*  o  to  t  and  hence  b(t) 
would  have  a  positive  Halting  value* 


We  note  that  in  Theorem  2,  since  F(x,t)  is  assumed  to  be  a  proper 
probability  distribution  function*  it  is  trivially  true  that  lia  f (x,t)  ■ 

11m  •=■?  F(x,t)  -  o. 

t-*»  3t 

6.  Susaary 

In  this  report  we  have  established  a  splitting  theoren  for  NHPP  similar 
to  the  splitting  theorem  for  stationery  Poisson  Processes.  Using  this 
powerful  result  we  showed  that  the  output  process  of  an  M(t)/6(t)/«*  queue 
is  a  NHPP  and  that  the  lumber  of  departures  from  the  queue  during  a  finite 
Interval  are  Independent  of  the  number  of  customers  remaining  in  service  at 
the  end  of  the  Interval.  Furthermore,  using  the  splitting  theorem  we  showed 
that  the  number  of  busy  servers  in  the  M(t)/G(t)/«  queue  is  Poisson  distributed 
with  mean  /  A(x)F(x,t)dx.  These  results  provide  a  tool  for  analysing  queues 
for  which  She  arrival  process  and  service  times  are  non-s tat ionary. 


Hlllestad,  R.  J.  and  M.  J.  Carlllo,  "Mod* la  and  Tachniqu—  for  Recoverable 
Iten  Stoclcage  Whan  Danand  and  tha  Rap air  Procaaa  ara  Mon-Stationary  — 
Part  I:  Parfonaanee  Measures , "  RAMD  Kota  N-148Z-AP,  Kay  1980. 

Klalnrock,  L.t  "Quaualng  Systems,  Vol  1:  Theory,"  John  Wiley  &  Sons, 

1975. 

Kitchen.  J.  V.,  "Calculus  of  One  Variable,"  Addlaon-Wessely  Publishing 
Co.,  1968. 

MLrasol,  N.  M. ,  "The  Output  of  an  H/G /•  Quaualng  System  la  Polsaon, " 

Latter  to  tha  Editor,  Operations  Research  (11:2)  pp282-284. 

Par sen,  E.,  "Stochastic  Processes,"  Bolden-Day,  Inc.,  1962. 

loss,  S.  M.,  "Applied  Probability  Models  with  Optimisation  Applications," 
Holdan-Day,  Inc.,  1970. 

Wlddar,  D.  V.,  "Tha  Laplace  Trans  fern,"  Princeton  University  Press,  19A6. 


DISTRIBUTION 


Deputy  Under  Sec’y  of  the  Arm/,  ATTN:  Office  of  Op  Beech 
Headquarters,  US  Anr  Materiel  Development  &  Beedlneee  Coemend 
DRCDM-S 

Commandant,  US  Army  Logistics  Management  Center,  Ft.  Lee,  FA  23801 
Defense  Logistics  Studies  Info  Exchange,  DBXMC-D,  ALMC,  Ft.  Lee,  FA  23801 
Defense  Technical  Info  Center,  Cemeron  Stetlon,  Alexendrle,  FA  22314 
Commander,  USA  Armament  Materiel  Readiness  Cmd,  Bock  Island,  IL  61299 
DB5AR-SA 

Commander,  USA  Communlc atlona-Electronica  Cmd,  Ft.  Monmouth,  NJ  07703 
DRSEL-PL-SA 

Commander,  USA  Missile  Command,  Redstone  Arsenel,  AL  35898 
DRSMI-S 


Commander,  USA  Troop  Support  6  Aviation  Materiel  Reediness  Command, 

S».  Louis,  MO  63120 
DRSTS-BA 

Commander,  US  Army  Teak-Automotive  Cwanrt,  Warren,  MI  48090 
ATTN:  DRSTA-8 

Commander,  USA  Armament  RAD  ted,  ATTN:  DXDAR-SE,  Dover,  NJ  07801 
Cownder,  USA  Aviation  SAD  ted,  4300  Goodf allow  Blvd,  St.  Louis,  MO  6312 
Commander,  USA  Aviation  RAD  ted,  ATTN:  DXDAF-BD  (Marl  Glenn),  4300 
Goodf allow  Blvd,  St.  Louis,  MO  63120 
Commander,  USA  Electronics  RAD  Cud,  ATTN:  DRDKL-ST-SA,  Bn  2D315, 

Bldg.  2700,  Ft.  Monouth,  NJ  07703 

Commander,  USA  Mobility  Equipment  RAD  ted,  ATTN:  DRDME-O,  Ft.  Bel  voir, 

FA  22060 

Commander,  USA  Natick  RAD  ted,  ATTN:  DBXHK-O,  Natick,  MA  01760 
Director,  Army  Materiel  Systems  Analysis  Activity,  AFG,  MD  21005 
DRZST-0 


DBXST-R 


DRZST-FO 

DRXST-DA 

Commander,  US  Army  Logistics  Center,  Ft.  Lee,  FA  23801 
Commander,  US  Air  Force  Logistics  Cad,  NFAFB,  ATTN:  AFLC/XRS,  Dayton, 
Ohio  45433 

US  Navy  Fleet  Materiel  Support  Office,  Naval  Support  Depot, 
MechaalcOburg,  FA  17055 

Ooorga  Washington  Uaivarsity,  Inst,  of  Management  Science  A  Bugr . , 

707  22nd  St.,  N.V.,  Wash.,  DC  20006 

Naval  Foatgraduete  School,  ATTN:  Dept  of  Opns  teal,  Monterey,  CA  93940 


1 


1 


1 


1 


Air  force  Institute  of  Technology,  ATTN:  SLGQ  Heed  Quant ltat ire  Studies 
Dept.,  Dayton,  OH  43433 
US  Any  Military  Acadeay,  Heat  Point,  NT  10996 

Scientific  Advisor ,  ATCL-SCA,  Any  Logistics  Center,  Ft.  Lee,  VA  23801 
Librarian,  Logistics  Mgt  Inst.,  4701  Sangaaore  Rd. ,  Wash., DC  20016 
University  of  Florida,  ATTN:  Dept  of  Industrial  Systems  Engr. , 
Gainesville,  FL  32601 

HAND  Corp. ,  ATTN:  S.  M.  Dresner,  1700  Main  St.,  Santa  Monica,  CA  90406 
Director,  Any  Menageasnt  Engineering  Training  Agency,  Hock  Island 
Arsenal,  Rock  Island,  IL  61299 

Defense  Logistics  Agcy,  ATTN:  DLA-LO,  Cameron  Sta,  Alexandria,  VA  22314 
Dap  Chf  of  Staff  (I&L),  HQ  USMC-LMP-2,  ATTN:  LTC  Sosneborn,  Jr., 

Wash.,  DC  20380 

Operations  Research  Center,  3115  Etcheverry  Hell,  University  of 
California,  Berkeley,  CA  94720 

Dr.  Jack  Mucks t ad t.  Dept  of  Industrial  Engineering  &  Operations  Research, 
Upson  Hall,  Cornell  University,  Ithaca,  NT  14890 
Dept  of  Industrial  Engineering,  University  of  Michigan,  Ann  Arbor, 

MI  48104 

Prof  Robert  M.  Stark,  Dept  of  Stat  &  Computer  Sciences, 

University  of  Delaware,  Newark,  DE  19711 
Prof  E.  Gerald  Hurst,  Jr.,  Dept  of  Decision  Science,  The  Wharton 
School,  University  of  Penna.,  Phila.,  PA  19104 
Dept  of  Industrial  Engr.  &  Engr.  Management,  Stanford  University, 

Stanford,  CA  94305 

Commander,  US  Any  C omnications  Command,  ATTN:  Dr.  lorry, 

CC-LOG-LEO,  Ft.  Huschuca,  AZ  85613 
Commander,  US  Army  Test  6  Evaluation  Cad,  ATTN:  DRSTE-ST, 

Aberdeen  Proving  Ground,  MD  21005 
Prof  Harvey  M.  Wagner,  Dean,  School  of  Business  Ada,  University 
of  North  Carolina,  Chapel  Hill,  NC  27514 
DARG0M  Intern  Training  Center,  Rad  River  Any  Depot,  Texarkana,  TX  75501 
Prof  Leroy  B.  Schwarz,  Dept  of  Management,  Purdue  University, 

Krannert  Bldg,  West  Lafayette,  Indiana  47907 
US  Any  Training  A  Doctrine  Comuad,  Ft.  Monroe,  VA  23651 
Operations  6  Inventory  Analysis  Office,  NAVSUP  (Code  04A)  Dept 
of  Navy,  Wash.,  DC  20376 

US  Any  Research  Office,  ATTN:  Robert  Launer,  Math.  Div., 

P.0.  Box  12211,  Research  Triangle  Park,  NC  27709 
Prof  William  P.  Plerakalla,  3641  Locust  Walk  CE,  Phila.,  PA  19104 
Air  Force  Logistics  Mgt  Center,  Gunter  Air  Force  Stetlon,  AL  36144 
AFLMC/LGT 

AFLMC/XRP,  Bldg.  205 

Engineer  Studies  Center,  6500  Brooks  Lana,  Wash.,  DC  20315 
Dr.  Carl  Walsmsn,  C.A.C.I.,  1815  N.  Fort  Igrer  Drive,  Arlington,  VA  22209 
Commander,  US  Any  Missile  CM,  ATTN:  Ray  Dotson,  DRSKI-DS,  Redstone 
Arsenal,  AL  35898 

Prof  Peter  Robot,  Graduate  School  of  NpU,  University  of  Rochester, 
Rochester,  NT  14627 

Prof  Barney  Blsslnger,  The  Pennsylvania  State  University,  Middletown, 

PA  17057 

Prof  Garry  D.  Sc udder,  Mgmt  Sciences  Dept,  College  of  Business  Adm, 

University  of  Mlnneaote,  271  19th  Ave.,  South,  Mlnneepolls,  MN  55455 


18 


